% 频率采样法
close all;
clear;
N=30;
H=[ones(1,4),zeros(1,22),ones(1,4)];
k=0:(N/2-1);k1=(N/2+1):(N-1);k2=0;
A=[exp(-j*pi*k*(N-1)/N),exp(-j*pi*k2*(N-1)/N),exp(j*pi*(N-k1)*(N-1)/N)];
HK=H.*A;
hn=ifft(HK);
fs=15000;
[c,f3]=freqz(hn,1);
f3=f3/pi*fs/2;
figure(1)
stem(hn);
figure(2)
plot(f3,20*log10(abs(c)));
title('幅频响应特性特性');
xlabel('频率/HZ');ylabel('衰减/dB');
grid;
